Lieb-Robinson bounds and the speed of light from topological order 
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We apply the Lieb-Robinson bounds technique to find the maximum speed of interaction in a spin model 
with topological order whose low-energy effective theory describes light [see X.-G. Wen, Phys. Rev. B 68, 
115413 (2003)]. The maximum speed of interactions is found in two dimensions is bounded from above less 
than v2e times the speed of emerging light, giving a strong indication that light is indeed the maximum speed 
of interactions. This result does not rely on mean field theoretic methods. In higher spatial dimensions, the 
Lieb-Robinson speed is conjectured to increase linearly with the dimension itself. Implications for the horizon 
problem in cosmology are discussed. 
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Introduction. — The principle of locality is one of the most 
fundamental ideas of modern physics. It states that every 
physical system can be influenced only by those in its neigh- 
borhood. The concept of field is the outcome of taking this 
principle seriously: if object A causes a change on object B, 
there must be changes involving the points in between. The 
field is exactly what changes. In addition, if something is 
"happening" at all the intermediate points, then the interac- 
tion between the objects must propagate with a finite speed. 
Relativistic quantum mechanics is built by taking the locality 
principle as a central feature. In non-relativistic quantum me- 
chanics the situation is more subtle: signals can propagate at 
every speed and quantum correlations are non-local in then- 
nature. One can, in fact, send information over any finite dis- 
tance in an arbitrary small time [1]. However, the amount of 
information that can be sent decreases exponentially with the 
distance if the Hamiltonian of the system is the sum of local 
pieces. Specifically there is an effective light cone resulting 
from a finite maximum speed of the interactions in quantum 
systems. This is the essence of the Lieb-Robinson bounds 
This notion have recently attracted interest in the context of 
quantum information theory, condensed matter physics, and 
the creation of topological order HI [1,3,151]. 

The concept of topological order is one of the most pro- 
ductive recent ideas in condensed matter theory [6]. It pro- 
vides explanations for phases of matter (for example, frac- 
tional quantum Hall liquids) that cannot be described by the 
paradigm of local order parameters and symmetry breaking. 
If local order parameters cannot describe such phenomena, 
then their order could be of topological nature [6] . Topologi- 
cal order gives rise to a ground state degeneracy that depends 
on the topology of the system and is robust against any lo- 
cal perturbations [7]. Because of this property, topologically 
ordered systems appear to be good candidates for robust quan- 
tum memory and fault- tolerant quantum computation |8]. 



Not only can topological order explain exotic phases of 
matter, but it offers a whole new perspective to the problem 
of elementary particles. There are particles that we regard as 
fundamental, like photons and fermions, and other particles 
that can be interpreted as collective modes of a crystal. For 
example, we can describe phonons in this way because of the 
symmetry of the crystal. The understanding of the phases of 
matter provides an explanation for the phonon and other gap- 
less excitations. However, one can also ask whether photons, 
electrons, or gravitons are emergent phenomena too, not el- 
ementary particles. Let us consider the case of light. Pho- 
tons are U(l) gauge bosons and they cannot correspond to the 
breaking of any local symmetry [9[]. Nevertheless, they can 
be collective modes of a different kind of order, and this is 
the case of topological order. Indeed models with topologi- 
cal order can feature photons, fermions and even gravitons as 
emerging collective phenomena fioll . 

Light emerges from topological order as the effective low- 
energy theory of a quantum spin system. The quantum spin 
system is built as a local bosonic model, namely a system in 
which the principle of locality is enforced by the fact that the 
Hilbert space decomposes in a direct product of local Hilbert 
spaces and all the observables have to commute when far 
apart. Moreover, the Hamiltonian must be a sum of local 
observables. In the low-energy sector, and in the continuum 
limit, the effective theory can be described by the Lagrangian 
of electromagnetism. Therefore low-energy excitations be- 
have like photons. Maybe this is what photons really are, col- 
lective excitations of a spin system on a lattice with Planck- 
scale distance. But then, why do we not see signals that are 
faster than light? There could be all sorts of interactions that 
can propagate as fast as permitted by the coupling constants 
of the underlying spin model. A theory of light as an emer- 
gent phenomenon needs to explain why we do not see signals 
faster than light. 



In this letter, we exploit the Lieb-Robinson bounds to show 
that the maximum speed of the interactions is of the same or- 
der of magnitude than the speed c of emerging light. This an- 
swers why we can think of light as an emergent phenomenon 
and still not see any faster signals in this model. In the last 
part of the paper, we argue that the maximum speed increases 
linearly with the dimension of the space and consider the im- 
plications for the horizon problem in cosmology. 

Topological Order and Artificial Light. — If we want to im- 
pose the principle of locality in a strong sense, we must con- 
sider local bosonic models |90- Fermionic models are not re- 
ally local because fermionic operators do not generally com- 
mute even at distance. A local bosonic model is a theory 
where the total Hilbert space is the tensor product of local 
Hilbert spaces, local physical operators are finite products act- 
ing on nearby local Hilbert spaces, and the Hamiltonian is a 
sum of local physical operators. Thus local physical oper- 
ators must commute when they are far apart. If we restrict 
ourselves to the case of a discrete number of degrees of free- 
dom and finite-dimensional local Hilbert spaces, we have a 
quantum spin model. A quantum spin model can be therefore 
defined as follows. To every vertex X in a graph G we asso- 
ciate a finite dimensional Hilbert space Ti x . The total Hilbert 
space of the theory is H — <8xeG%£- To every finite subset 
of vertices X C G, we associate the local physical operators 
with support in X as the algebra B(Hx) of the bounded lin- 
ear operators over the Hilbert space Tix — ®igx^i- The 
Hamiltonian will have the form H\ oca i — J^xcG ®X> where 
to every finite subset X C G we associate an hermitian opera- 
tor $x with support in X. An example of local bosonic model 
is given by a spin 1/2 system on a lattice. To every vertex x 
in the lattice we associate a local Hilbert space fi x — C 2 . Lo- 
cal physical operators are finite tensor products of the Pauli 
matrices at every vertex. 

The bosonic model we consider is a lattice of quantum ro- 
tors. Its low-energy effective theory is a U(l) lattice gauge 
theory whose deconfined phase contains emergent light. Con- 
sider a square lattice whose vertices are labeled by i, with an- 
gular variable 6-^ and angular momentum S'S on its links. The 
Hamiltonian for the quantum rotor model is given by 



H rotor = c/e(E 5 ^) +jJ2^ s u 2 

i \ a / i,a 

+ E (t< aiQ2 >e l ( e ^- e ^ + h.c.) , 

!.{<*!. a 2 } 
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where a = ±1/2(1, 0), ±1/2(0, 1) are the vectors of length 
1/2 pointing towards the lattice axes [0]. In the limit t, J <C 
U, the first term of the Hamiltonian H rotor behaves like a 
local constraint and makes the model a local gauge theory. 
Defining g :— 2/U(t\2t-i-2 + ti-it-ii), the effective low- 




F1G. 1: (Color online) A 2D— dimensional rotor lattice. To every 
plaquette p is associated a rotor operator W v as a function of the 
variables 8ij. The graph G is the one drawn in thin black lines. The 
graph G' is the graph with black and blue (lighter, bigger) dots as 
vertices and blue thin lines as edges. The red dashed line shows a 
path of length n — 22 from the point P to the point Q which are at 
a distance 2d(P, Q) — 8 on G' or d(P, Q) = 4 on G. These paths 
contain alternating link and plaquette operators. 



energy theory becomes 

ffe// = -/E( 5 u) 2 -sE 
<ij> P 

where W n = e l(ei2 ~ e23+834 - 



W + h.c. 



E 

<ij>,P 



(1) 



p — c. ■ 41 ' is the operator that creates a 

string around the plaquette p (see Fig. [D and the t's are cou- 
pling constants. Although a lattice gauge theory is not a local 
bosonic model, this does not violate locality because H e f f is 
just an effective theory. The fundamental theory is local and 
H e f f is still a sum of local terms. In the large g/J limit, the 
continuum theory for the Hamiltonian H e tf is the Lagrangian 
of electromagnetism 
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with speed of light given by c = \J2gJ. 

Lieb-Robinson Bounds and the speed of sound in spin 
systems. — Here we review the proof of the standard Lieb- 
Robinson bounds J2|] in the variant first proven in 04J and 
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also exposed in J3l. We consider a Hamiltonian of the form 
H := J^xcG Now consider an operator Oy with sup- 
port in a set Y C G. The time evolution for this operator 
under the unitary induced by H is Oy(t) = e ltH Oye~ ltH . 

The Lieb-Robinson bound is an estimate of an upper bound 
of the commutator of two operators Op(t), Og(i') with sup- 
port in different regions P and Q and at different times t 
and t'. If the interaction map <£>x couples only nearest- 
neighbor degrees of freedom, the Hamiltonian can be written 
as H = X)<ij> an d tne Lieb-Robinson bound reads 

\\[O P (t),O Q (0)]\\ < nO P \\\\0 Q \\f^ mk ™* r N PQ (n) 

< 2\\Op\\\\0 Q \\Cexp[-a(d(P,Q)-vt)} 

where h max — max < ij >6( 3 /iy and Npg(s) is the number 
of paths of length s /2 between the points P, Q at distance 
d(P, Q) in G [1]. The constants C, a, v have to be determined 
in order to get the tightest possible bound. This bound is loose 
for several reasons: the crude maximization over fty, the over- 
look about the Hamiltonian's details, and the fact that all in- 
teractions are summed in modulus instead than amplitude, so 
that destructive interference is not taken in account. 

Lieb-Robinson Bound for the emergent U(X) model. — 
What do the Lieb-Robinson bounds tell us about the model 
H e f f with emergent light? Is the maximum speed of the in- 
teractions something like the speed of the emergent light or 
something completely different? As we have seen, this is of 
great importance if we want to take seriously the theory of 
light as an emergent phenomenon. 

If we apply naively the Lieb-Robinson bounds to the 
Hamiltonian of the U(l) lattice gauge theory, we see that the 
speed v is proportional to the strongest of the coupling con- 
stants, v (x g. Since light only exists in the phase g 3> J, 
we would have v ^> \fg~J . Fortunately, the bound can be 
made much tighter by examining the details of the Hamilto- 
nian and the specific way the interactions propagate. Consider 
the function /(i) := [Op(t), Oq(0)]. Then consider the set 
Zi := {Z C G : [$z,O p ] ^ 0}, the support of the com- 
plement of the commutant of Op in the set of interactions. It 
turns out [3] that /(<) obeys the differential equation 

/'(*) = -< E (m),* z (t)}+iop(t),[$ z (t),o Q (o)]}), 

ZCZ! 

where $z{t) = e $ze~ l . From the above equation, and 
using the norm-preserving property of unitary evolutions, we 
can establish [3] the bound 

I|[Op(*),O q (0)]||<||[O p ,Oq]|| 
+2\\0 P \\ f || £ [<fz(t),O Q (0)}\\. 

Successive iterations of the above formula yield 

\\[Op(t),O Q (0)]\\ <2||Op||||0 Q ||E i ^r a "' (2) 

n=0 



where 

n 

YiCZi z=l 

and we define Z i+1 := {Z C G : [$z, $z'czj ^ 0} and 
where 0\ and O2 are two non-commuting local operators of 
the Hamiltonian. The meaning of the above expression is the 
following. Every element of the sum is a product of the type 
Y\i ||$i|| such that [<5>i, 7^ for every i. If each <J>i is a 

local bosonic operator, every one of those products is a path 
on the lattice. 

Let us apply these considerations to the case of the effec- 
tive Hamiltonian H e ff- For the sake of simplicity, consider 
Op, Oq to be the spin operator S z at the points P, Q. For this 
Hamiltonian, the only non commuting operators are W p and 
Sf when they have a vertex in common (see FigJTJ. There- 
fore, a path in (O will consist of steps from a plaquette to any 
of the four links bordering it, alternated with steps from a link 
to any of the two incident plaquettes. Any such path is then a 
path drawn with dashed edges in Figfl]on the lattice G 1 . To 
every path of length n on G' will then correspond an operator 
whose norm is 

nr=ill $ ill = (g J ) n/2 - Therefore, denoting 
by N' P q(ti, d) the number of paths of length n on G' from P 
to a given point Q at a distance 2d, we obtain the following 
bound 

a n <N' PQ (n,d)(gJ)% (4) 

A gross bound is given by N' PQ {n,d) < 2V&"' e K(n ~ 2d+4: \ 
for every n > 0. This is because there are 8 ways to do a suc- 
cession of two steps on G'\ 4 choices from a blue vertex and 2 
from a black one (see FigJT), and this quantity is greater than 
one iff there is at least one path of length n between P and Q. 
Moreover, the iteration of Eq.© can be built by replacing the 
2|t| with Ip 1 ,* 2 ]!!^* 1 !!!!^ 2 !!)- 1 ^ = s/2\t\, and obtain 

\\[Op{t),0 Q {<m\ < 4e 4K ||0 P ||||0 Q || e - 2K ( d -^ Z£!1 l t l). 

Optimizing for k we get v^p = ^/2ey/2gJ = sple x c. 
We verify numerically that our approximation for N'(n,d) 
is good: An exact combinatorial formula is (we drop the sub- 
script PQ from now on): 

We numeric ally studied the quantity X) ^Lo ( ^\ ^ I )" On ! be- 
cause that is the one that enters the bound Eq. (f2j). The facto- 
rial at the denominator makes the series converge rapidly and 
we obtain, together with Eq. (f2]i 

\\[Op(t),O Q (0)}\\<nOp\\\\O Q \\Ae-^ (5) 

The speed v is estimated numerically as v w s/2e^/2gj = 
vlr — \/2e x c. Let us try to understand this result. Eq. 
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(0 establishes that all the observables that are outside of the 
effective light cone centered on P with speed of light vlr will 
have an exponentially small commutator with the observables 
in P. This result sets a limit to the speed of interactions in the 
spin system. It proves that any signal outside of a light cone 
generated with a speed that is of the same order of magnitude 
(and with the same dependence on coupling constants) of light 
will be exponentially suppressed. We consider this result a 
strong indication that the maximum speed of signals is light. 
So the theory of emerging light explains why its speed is also 
the maximum speed for any signal at low energies. If we were 
able to probe energies of order U we could still find faster 
signals. 

The cosmological horizon problem. — The isotropy of the 
cosmic microwave background presents us with the horizon 
problem: how is it possible that regions that were never 
causally connected have the same temperature? The hori- 
zon problem arises from the stipulation that interactions can- 
not travel faster than a finite speed, which defines a causal 
cone. Inflation solves the horizon problem by introducing 
an exponentially fast early expansion which allows for initial 
causal contact and thermalization of the observable universe 
ifTHl . Alternative proposed solutions require a mechanism for 
changing the speed of light as we trace the history of the uni- 
verse backwards in time [12] or a bimetric theory 111 311 . Dy- 
namically emerging light could also resolve the horizon prob- 
lem. We now wish argue that our results on the maximum 
speed of interactions in speed systems can be used to justify 
this. 

Let us first understand how the speed vlr depends on 
the dimension D of the space. Consider a hypercubic 
D— dimensional lattice. The number of paths of length n 
on a hypercubic D— dimensional lattice will be Nr>{n) ~ 
[AD(D - l)]"/ 2 . If the two dimensional case is any indica- 
tion, a good enough approximation for Nrj(n, d), the number 
of paths of length n between two points P and Q at distance 
2d apart will thus be N D (n, d) ~ [AD{D - i)]»/V ( "~ 2d) , 



which implies a speed vlr{D) growing linearly with D IU5I1 . 
Now consider a model of the universe in which we start with 
an extremely connected graph that evolves towards a less and 
less connected graph, for instance a hypercubic lattice of di- 
mension D{t) = Di n (l — at). This type of situation has 
been hypothesized in quantum spin models of the universe 
like quantum graphity II 1411 . In such a system, the maximum 
speed of interactions will decrease linearly in with the dimen- 
sion D of the space, and hence in time, providing a possi- 
ble explanation to the horizon problem in cosmology. Light 
cones then have parabolic sides and allow correlations at early 
times without violating causality since the distance of corre- 
lated points is of the order of a(U — tf) 2 . 

Conclusions. — In this letter, we applied the technique of 
the Lieb-Robinson bounds to estimate the maximum speed of 
interactions for a quantum spin model with topological order 
and emerging U(l) gauge symmetry. The importance of this 
model is that the low energy excitations are photons. Light can 



be regarded as an emergent phenomenon and photons can be 
seen as collective modes instead of elementary particles [9]. 
This theory poses the problem of why we do not see other ex- 
citations that are faster than light. The technique of the Lieb- 
Robinson bounds, in the variation presented here, shows that 
the maximum speed of excitations in the model has the same 
order of magnitude as the speed of light. Of course, it is easy 
to construct a different model where there is emerging light 
and other faster particles. One of the fundamental questions 
of physics is to explain why this does not seem to happen in 
nature. In order to address this question, the Lieb-Robinson 
technique could prove useful if it could be modified in order to 
find a tight bound for a frustrated model, where destructive in- 
terference prohibits other signals potentially faster than light. 
In perspective, we think that this technique can prove useful to 
find exact results in 2D condensed matter models where there 
is scarcity of results that are not just numerical. 

Finally, we have discussed the implications of the finite 
speed of signals for cosmology and the horizon problem. 
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